Annihilations of superheavy dark matter in superdense clumps 
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Superheavy dark matter (SHDM) exchanges energy with its environment much slower than parti- 
cles with masses close to the electroweak (EW) scale and has therefore different small-scale clustering 
properties. Using the neutralino as candidate for the SHDM, we find that free-streaming allows the 
formation of DM clumps of all masses down to ~ 260 m x in the case of bino. If small-scale clumps 
evolve from a non-standard, spiky spectrum of perturbations, DM clumps may form during the ra- 
diation dominated era. These clumps are not destroyed by tidal interactions and can be extremely 
dense. In the case of a bino, a "gravithermal catastrophe" can develop in the central part of the 
most dense clumps, increasing further the central density and thus the annihilation signal. In the 
case of a higgsino, the annihilation signal is enhanced by the Sommerfeld effect. As a result anni- 
hilations of superheavy neutralinos in dense clumps may lead to observable fluxes of annihilation 
products in the form of ultrahigh energy particles, for both cases, higgsinos and binos, as lightest 
supersymmetric particles. 
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PACS numbers: 12.60.Jv, 95.35.+d, 14.80.Nb, 98.70.Vc, 95.85.Pv 



I. INTRODUCTION 

The case for the existence of non-relativistic, non- 
baryonic dark matter (DM) in the universe is stronger 
than ever [lj. But although a wealth of observational 
data provides compelling evidence for a ~ 23% contribu- 
tion of DM to the total energy density of the universe, its 
nature is still not known. The most popular DM type are 
thermal relics, i.e. particles that were at least once during 
the history of the Universe in chemical equilibrium with 
the thermal plasma. 

The present relic abundance il x of a thermal relic 
scales approximately with its annihilation cross section 
Moreover, unitarity bounds anni- 
hilations as a 
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J 2 and thus the observed value [2E] 



f^CDM^ 2 = 0.1 of the DM abundance constrains the an- 
nihilation cross section as (cr ann u) ~3x 10 _26 cm 3 /s and 
limits the mass of any thermal relic as m x < 50 TeV [|| . 
Hence thermal relics offer detection prospects both for 
direct and indirect searches as well as at accelerator ex- 
periments. 

The assumption that the DM particle was once in 
chemical equilibrium is however not necessary and does 
not hold in particular for sufficiently heavy particles. Su- 
perheavy particles are generated at the end of inflation 
and they can play the role of DM particles 0, Q ■ Grav- 
itational production at the end of inflation provides the 
most natural mechanism for the generation of superheavy 
dark matter (SHDM) @. Their decays can result in the 
observable signal in the form of UHE gamma-rays [ij and 
UHE neutrinos [7], if they are metastable and long-lived. 

While the idea of SHDM is theoretically appealing, 
the detection of SHDM is challenging. Clearly, accel- 
erator searches and direct detection are not feasible in 
the case of SHDM. The feasibility of indirect detec- 



tion of stable SHDM depends on its annihilation rate 
N ann oc (p/m x ) 2 (<T ann u) that in turn scales naively as 
•^Vann °c m x^' Since backgrounds like cosmic rays from 
astrophysical sources or the diffuse photon flux decrease 
only as 1/E a with a < 3, indirect detection of DM seems 
to become more and more difficult for increasing DM 
masses. The only possibility which overcomes this diffi- 
culty is the superdense central region of DM clumps [|| , 
but one needs the realistic scenario for the very high den- 
sity of DM in the clump center or formation of superdense 
clumps UBS. 

We shall use as candidate for SHDM the neutralino in 
the model of superheavy supersymmetry, as suggested in 
Ref. [ll| . Superheavy supersymmetry is a unique scheme 
that respects perturbative unitarity despite of coupling 
particles with mass much larger than the weak scale to 
the electroweak sector. Within this model, the lightest 
supersymmetric particle (LSP) that we choose as the neu- 
tralino is a natural candidate for SHDM. 

Aim of the present work is to study the detection 
prospects for stable SHDM. Since the annihilation signal 
from the mean distribution of DM in the halo is far below 
observational limits, we examine if new effects specific for 
SHDM exist that can improve the detection chances. 

One such effect can result from the early kinetic decou- 
pling of SHDM: While the mass spectrum of DM clumps 
formed by standard neutralinos with mass in the 100 GeV 
range has a cutoff at M min - (10~ 12 - 10- 4 )M Q Hid, 
the cutoff can be diminished significantly e.g. in the case 
of ultra-cold WIMPs [10, 0- We will show that the 
cutoff is practically absent for a superheavy bino as DM 
particle, and clumps of all masses are possible beginning 
formally from ~ 260 m x . This low-mass cutoff increases 
the diffuse flux of UHE particles produced by annihila- 
tions. 
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Another effect is the formation of superdense clumps, 
in which the annihilation rate can be strongly enhanced. 
The formation of superdense clumps is discussed in gen- 
eral in the accompanying paper I (l5| . Here we study 
this problem in more detail for superheavy particles. 

The initial mass spectrum of DM clumps is determined 
by the spectrum of cosmological density perturbations. 
The simplest models of inflation predict a nearly scale- 
invariant power-law form that is then normalized to the 
COBE observations. Clumps formed in this case are not 
very dense, and the small-scale structure of DM enhances 
the annihilation signal not strongly. A quite different 
scenario arises, if the potential V(<f>, T) of the inflaton 
field contains features like a zero derivative V'(<j),T) for 
some field value <fi. The spectrum of perturbations in this 
case has a spike and clumps of SHDM can form already in 
the radiation-dominated (RD) epoch, and they can have 
very large densities. 

For sufficiently dense clumps, the relaxation processes 
due to gravitational two-body scatterings can initiate a 
"gravithermal catastrophe" and the density profile of the 
clump core steepens to an isothermal profile p oc r~ 2 
with a much smaller new core radius. While this process 
was discussed in general in Ref. (hereafter Paper I), 
we show in this work that for a bino this process can 
take place and determine the required initial conditions. 
Annihilations of DM in such dense clumps are strongly 
amplified because of the density enhancement. 

In contrast to a bino, winos and higgsinos are stronger 
coupled to the thermal plasma. As a result, a "gravither- 
mal catastrophe" does not develop. However, the veloci- 
ties of DM particles in superdense clumps are very small, 
leading to an enhancement by the Sommerfeld effect in 
case of higgsinos or winos. Taking into account both ef- 
fects, we find that annihilations of SHDM bound in such 
superdense clumps may lead to observable fluxes of ul- 
trahigh energy particles for all types of neutralinos. 

Note that (in contrast to the case with standard power- 
law spectrum of cosmological perturbations) superdense 
clumps produced in the RD epoch from a spiky spectrum 
are not destroyed by tidal forces and their mass function 
peaks near a definite value. Therefore the fraction of DM 
in the form of such clumps is £ ~ 1/2. 

This article is organized as follows. We start with a 
discussion of the energy relaxation time of SHLSPs in 
Sec. [TTl where we also derive the minimal mass of the 
SHDM clumps considering free-streaming and the hori- 
zon scale at kinetic decoupling. In Sec. IHII we summarize 
the evolution of the density profile of superdense clumps. 
Next we discuss the prospects to detect SHDM clumps 
in Sec. II VI and conclude in Sec. [V] 

We shall use below the following abbreviations: SHDM 
for superheavy dark matter and SHLSP for superheavy 
lightest supersymmetric which in our case is the super- 



heavy neutralino We fix the mass m x of the SHLSP 



as m x — 10 GeV, i.e. within the range suggested by 
gravitational production with fl x h 2 ~ 0.1 at the end of 
inflation. 



II. ENERGY RELAXATION AND THE 
MINIMAL MASS OF DM CLUMPS 

A. Relaxation time 

We shall consider here the energy relaxation of SHDM 
particles interacting with a thermal background at tem- 
perature T. The consideration is relevant for both ordi- 
nary and superdense clumps. 

The energy exchange between su perh eavy neutralinos 
and light fermions was calculated in |ll| for temperatures 
T below the weak scale, i.e. in the limit T <C mz <C 
Msusy with Msusy = min{Mi, M2, /i}. Here, /i de- 
notes the Higgs mixing parameter, and Mi and M2 the 
U(l) and SU(2) soft SUSY breaking masses. At lowest 
order in m 2 z / Mj USY , the neutralino masses are simply 



{Mi,M 2 , — for more details see e.g. [11 1 



The energy relaxation time is calculated as 
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2E k m x ^ J J 



... / da e l,i 



(1) 

where E k — (3/2)T is the mean kinetic energy of the 
neutralinos, Sp the neutralino momentum obtained in one 
scattering, 



(Sp) 2 = 2u; 2 [l-cos(6<)] 



(2) 



and the number density of relativistic fermions or bosons 
with g polarization degrees and energy uj is 



9, 



2tt 2 e"l T ± 1 



2tt 2 



uj 2 e-"' T . 



(3) 



Kinetic decoupling occurs when the energy relaxation 
rate t~T becomes smaller than the expansion rate H 
of the universe. During the radiation dominated (RD) 
epoch, H = l/(2i), and 



H = 1.66^5* 



Mpi ' 



(4) 



where denotes the number of relativistic degrees of 
freedom and Mp\ the Planck mass, Mpi = 1/\/Gn — 
1.2 x 10 19 GeV. 

As the calculations of Ref. [ll| show the energy- 
relaxation time due to elastic scattering of neutralino on 
the background fermions at T < mz is larger than the 
corresponding age of the universe. Thus, it is necessary 
to include also the elastic scattering of neutralino on weak 
gauge bosons and the light higgs for mz «T« Msusy- 
Since the mass splitting between the lightest neutralino 
and chargino is of order 0(m 2 z /Msusy) for a higgsino or 
wino, inelastic processes like x° + v c — > X + + e ~ con- 
tribute also to the energy exchange between neutralinos 
and the plasma. As we will see, elastic interactions with 
fermions are for T 3> mz subdominant, while inelastic 
ones give the dominant contribution to the energy relax- 
ation of the neutralino. Since the latter are absent for 
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a bino, the relaxation time is strongly dependent on the 
nature of the neutralino. We consider in the following 
only the cases of a bino and a higgsino, noting that a 
wino behaves in most respects as a higgsino. 



1. The Bino as the LSP 

The processes dominating the energy exchange be- 
tween superheavy binos and the thermal plasma at T 3> 
mz are x% — > X%> %M /± — > xW^ and x^ 1 xh- In 
the first process only light higgs exchange contributes at 
leading order and the squared amplitude is 



\M x z- 



2 _ e A M^{n%m{2P) + M 1 ) 2 
' xZl ~ 3cos 4 (0vtO(M 2 --^i) 2 



(5) 



where tan (3 — v% jv-i is the ratio of the two Higgs vev 
and Qw the Weinberg angle. The matrix elements of 
the other relevant processes have the same form, and the 
total matrix element squared weighted by the relevant 
polarization degrees is given by \M.\ 2 = ^Zi9i\M-i\ 2 — 
l2\M x z^ x z\ 2 ■ The resulting energy relaxation time fol- 
lows as 



Trcl 



32a 2 T 4 (Af 1 +/is 2(3 ) 2 



(6) 



In the last step, we simplified the energy relaxation time 
T rc i introducing a common mass scale MgusY and A = 



TIC 



/(32a 2 



w 



2. The Higgsino as the LSP 

The inelastic process x\ v & ~> Xi e ~ 1S dominated by W 
exchange. Its leading contribution to the squared matrix 
element is given by 

= 2e Videos 2 (0/2) 
' x**-*^-' ^(^(I-oob^J+O 3 ' 

The energy relaxation time becomes for /i> T> mw 

67rsin (6w) A* 



Trol 



a 2 7V eff T 2 [In (4T 2 /m 2 w ) - 2 7 + 1] 



and 



Trcl 



AWT 6 



with B 



160a 2 



(8) 



(9) 



for m x ± — m x o <C T <C mw The factor N e g counts the 
number of fermions contributing to the inelastic processes 
at temperature T. 



B. Minimum mass of DM clumps 

The growth of small-scale fluctuations can be smeared 
out by destructive processes such as free-streaming, 
acoustic oscillations and others [H, [l6| (see appendix A 
in Ref. [2(| for a general discussion). For neutralinos 
with mass close to the electroweak (EW) scale, these 
processes determine the minimal wave-length in the per- 
turbation spectrum that can grow and thereby also the 
minimum clump mass. Here we consider SHLSPs and 
study the effectiveness of these damping effects in the 
case m x ^> mz- 

Let us calculate the cosmological age td and the tem- 
perature Td of kinetic decoupling of SHLSPs from the 
cosmic plasma, i.e. the moment td when the relaxation 
rate t~, equals the expansion rate H(td) of the Universe. 
In evaluating the condition for decoupling, ~ H, 
we set N c g — 431/4 = 107.75 as number of relativis- 
ts degrees in the standard model (SM) for T > m t 
and use for the running of coupling and mixing parame- 
ters with temperature T the SM relations, sin 2 9w{T) = 
l/6 + 5a(T)/[9a s (T))]. For M SU SY = 10 12 GeV we ob- 
tain then 



2 x 10 n GeV, 
2 GeV, 



bino 
higgsino 



(10) 



as decoupling temperature for a bino and higgsino, re- 
spectively. In the former case, <?* = 298/4 and N g = 66. 
For t > td, the SHLSP is not longer in thermal equilib- 
rium with the cosmic plasma and its momentum scales 
as p oc 1/a 2 . 

We consider now the physical processes relevant for a 
spherical region containing DM with total mass M close 
to the time of horizon crossing. The mass of DM inside 
the horizon as function of the temperature is given by 

M = 3.4 x 10 16 (T/100GcV)- 3 (iV cff /100)- 3/4 g. (11) 

In particular, at the temperature of kinetic decoupling 
given by Eq. (|10p the corresponding mass is equal to 



M d 



6 x 10~ 12 g, 
6 x 10 21 g, 



bino 
higgsino 



(12) 



For a bino, the mass Md is only 34 times greater than 
the particle mass m x - 10 11 GeV= 1.78 x 10~ 13 g. 

The evolution of fluctuations with mass M <C Md and 
M ^> Md is rather different after horizon crossing [l6| . 
Fluctuations in DM with mass M <C Md run out as sound 
waves in the radiation plasma. These fluctuations do not 
have a kick in the peculiar velocities of their DM particles 
and therefore do not grow logarithmically. After kinetic 
decoupling their amplitude freezes in until the matter 
dominated (MD) epoch, and their evolution is analogous 
to the evolution of entropy perturbations and described 
by the Meszaros solution (see e.g. 0). Therefore, there 
is a steepening in the mass spectrum below M ~ Md- 

In the opposite case, M ^> Md, the peculiar velocities 
just after the horizon crossing are equal to v p h — Shc/3 
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17). In contrast to thermal velocities these peculiar ve- 
locities are regular and directed toward the center of the 
fluctuation. The fluctuations grow according to the adi- 
abatic law S cx ln(t) due to the evolution of the peculiar 
velocities as v p (t) ~ v p #a(i#)/a(i). The free streaming 
scale Af s of SHDM is very small. Expressed in comoving 
units, this scale is growing during the RD epoch as 

\f s {t) = a(t ) / — — — = 2t d v d In , (13) 

Jt d a \t) a d a-d 

where v(t) = Vda(td)/a(t), Vd = (STd/m) 1 / 2 , and the 
ratio of the scale factors a(t)/a{td) is calculated from the 
Friedmann equation. The corresponding free streaming 
mass, 

47T 

M is (t) = —p c (t )n m , \i(t) , (i4) 

stops growing near the epoch of matter-radiation equal- 
ity, t ~ t cq . For the case of a bino in Eq. ([TU| . the 
decoupling time is td — 7 x 10~ 30 s, and we find 

1/4 1/4,3/2 / rri \ 3/2 r O/l ~1 
fS G 3/4 \ m J UG Peq t d j 

~ 4.6 x 1CT 11 g. (15) 

This value is only 260 times larger the particle mass. 
Thus the free-streaming mass of SHDM defines the cutoff 
in the mass spectrum. Formally, all clump masses are 
possible from M > 260m x . 

In the case of a higgsino, M{ s <C m x , and free- 
streaming plays no role for the evolutions of density per- 
turbations. 

Therefore, the two mass scales Md and Mf s may play 
the role of M m - ln . In the case of a bino, M{ s > Md and the 
steepening of the mass spectrum starts at M m j n ~ Mf s . 
In the case of a higgsino, Mf s is very small and M m i„ ~ 
M d . 



III. CLUMP STRUCTURE FOR STANDARD 
AND SPIKY PERTURBATIONS 

Let us consider first the formation and evolution of 
clumps of SHDM for a standard power-law spectrum of 
fluctuations. Clumps of SHDM are produced and evolve 
accordingto the usual hierarchical model, as described 
in Ref. [21| and Paper I, with the essential difference 
that the minimum clump mass is now the one derived in 
Sec. IIIBI With an accuracy sufficient for this schematic 
consideration, we can use M m i n of order of 260m x for the 
case of a bino. 

The basic features of this scenario for SHDM clumps 
are the same as for ordinary DM clumps: Most clumps 
are destroyed by tidal interactions in hierarchical struc- 
tures and the surviving clumps could be further destruc- 
ted in the Milky Way (see Paper I) . SHDM clumps in this 
scenario have a rather small density and SHDM particles 



with their small annihilation cross section are unobserv- 
able through their annihilation products. They can be 
detected only gravitationally as discussed in Paper I. 

Let us come now to the case of a spiky perturbation 
spectrum, or to any other case with fluctuations growing 
in the RD epoch. In these cases superdense clumps can 
be produced. Since in the RD dominated epoch large- 
scale structures are absent, there is no tidal destruction 
of small clumps. These clumps evolve as isolated objects 
without any merging. The first stage of evolution, the 
ordinary gravitational contraction, proceeds in the stan- 
dard way leading to a p(r) cx r - density profile and a 
large core with radius R c ~ (0.01 — 0.1)i?, where R is the 
clump radius. At this stage the core size is restricted by 
tidal forces [2l[, by a decreasing mode of perturbations 
[23( | or by phase density constraints according to Liou- 
ville's theorem [24|. The "gravithermal instability" sets 
in, when the density of the core reaches a critical value 
that is determined mainly by the mass of the DM parti- 
cle, and leads to an isothermal density profile p(r) cx r~ 2 , 
with a new small core, determined by the properties of 
SHLSP. This so-called "gravithermal catastrophe" occurs 
under the influence of two-body gravitational scattering 
in full analogy of this process to the one in globular clus- 
ters. We will discuss this process now in some detail. 

A. Gravitational relaxation and evolution of the 
clumps 

The gravithermal instability in globular clusters sets 
in due to two-body gravitational relaxation. This pro- 
cess can be the dominant one for the superdense clumps 
from SHDM particles. The other relaxation channel, EW 
scattering of superheavy neutralinos loses the competi- 
tion, because its cross section is proportional to m~ 2 and 
the interaction has a short range. On the other hand, 
the very high clump density provides the relaxation time 
irci.gr to be shorter than the age of the universe to- 

The two-body gravitational relaxation time deter- 
mined as *-J jgr = (l/E)(dE/dt) in energy space can be 
taken from calculations for globular clusters (see [25[ and 
references therein) as 

1 v 3 
W ~ 4^G 2 m x nln(0.4A0' { ' 

where v ~ (GAf/i?) 1 / 2 is the virial velocity, and N is 
the total number of particles in the clump, N = M/m x . 
From Eq. ([16|) one observes indeed that the relaxation 
time is inversely proportional the mass squared m 2 of 
SHDM particle and the density n of the core. The loga- 
rithmic term ln(0.4./V) takes into account the long-range 
character of the gravitational interaction. 

The relaxation time shorter than the age to of the uni- 
verse leads to the "gravithermal catastrophe" , which re- 
sults in an isothermal density profile p(r) cx r~ 2 . In 
this regime the main process responsible for the evolu- 
tion of the clumps becomes the evaporation of particles 
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from the core, and the following calculations are in full 
analogy with the globular clusters case (25|. The escap- 
ing particles have approximately zero total energy, and 
therefore the energy of the core is approximately con- 
stant. The process of the core evolution can be described 
by the rate of evaporation and the virial theorem, which 
can be written as 



N c /N c = -a/t reU 



R c /R c = 2N C /N C - E c /E c 



(17) 

where a sa 7.4 • 10 3 for a Maxwellian initial velocity 
distribution. Joint integration with logarithmic accuracy 
[ln(0.4iV) = const in Eq. (fTf) ] and the condition E c = 
const gives the time evolution for the core mass M c and 
radius R c : 

M c (t) = m x N c (t) = M c>i (l - (t - U)/t e ) 2 '\ (18) 



R c {t)=R^{l-(t-t l )/t e f/\ 



(19) 



where t e = 2/(7a)i ro i,g r .i — 40i rc i ;gr ^ , and the subscript 
i marks the values at the initial moment U of the clump 
formation. The time t e is less than the age of the Uni- 
verse for clumps above the dotted line in Fig. Q] for 
m x = 10 11 GeV as mass of the DM particle. Thus for 
clumps above the dotted lines, relaxation results in the 
"gravithermal catastrophe" producing an isothermal pro- 
file p oc r~ 2 with a tiny new core. It diminishes until the 
central density becomes sufficiently large and new pro- 
cesses, such as annihilations or the pressure of a degen- 
erate Fermi gas, enter the game (see below). 



B. Formation of the new core 

When the large initial core loses its stability and con- 
tracts under the two-body gravitational forces, the gravi- 
tational relaxation time is increasing. It is caused by the 
increase of the core density as 



p(t) CX [1- {t-ti)/t e 



-10/7 



(20) 



which follows from Eqs. (1191) . In principle, this phe- 
nomenon can stop the contraction and stabilize the core. 

After the core collapsed, the singular profile p a r -2 
extends formally down to some small radius R c . For the 
density profile p(r) = p c (r / R c )~ 2 , the relative core radius 
x c = R c /R is given by x c = (p/Spc) 1 / 2 , where p and p c 
are the mean and the maximal density of the clump. Do 
any physical processes exist that prevent the extremely 
large densities, i.e. a very small radius, of the new core? 

The first candidate for such process is given by the 
EW elastic scattering of SHDM particles i.e. by sclfin- 
teraction. The relaxation time for this process can be 
estimated as 



4(2tt) 



1/2, 



rel,xx 



'XX' 



3 5 / 2 



(21) 



where a xx is the cross section of elastic neutralino scat- 



tering 



da 



xx 



Aa 2 



dn 



(22) 



where A is a constant of order one that depends on the 
SUSY parameters, a = 1/137, and v ~ (GM/R) 1 / 2 is 
the virial velocity. This relaxation time has the same de- 
pendence on the core-density (oc n c ) as the gravitational 
relaxation time (|16p. but is for SHLSPs many magni- 
tudes smaller. Therefore, self-interactions cannot stop 
the gravitational collapse. One may see this effect in a 
different way: The core remains transparent for super- 
heavy neutralinos down to extremely small radii. More- 
over, Q proved that elastic scatterings do not change 
the central distribution of DM. The effect of elastic relax- 
ations in the models of self-interacting DM was studied 
also, e.g., in [26| . 

Another effect which can stop the gravitational con- 
traction is SHDM particle annihilation. This effect was 



studied in Refs. [27J and [28j. In the former the core 



radius was found from the condition that the character- 
istic annihilation time in the core should be longer than 
the time of core formation estimated as hydrodynamical 
free-fall time ty L ~ [Gp)^ 1 ! 2 . For an isothermal density 
profile, the corresponding dimensionless core radius x c is 
given by 



1/2 



(23) 



The mass of superheavy particle m in denominator makes 
the core radius rather small, 



7.4-10 



-13 



IT) 



-3/2 
11 



10 5 g cm 



-3 



1/4 



(24) 



where m n = m/(10 n GeV). 

In a more reliable approach [28| , the radius of the core 
was estimated equating the core accretion rate with the 
annihilation rate inside it. The accretion rate was cal- 
culated from the Euler and Poisson equations. The cal- 
culated core radius is much smaller than the one from 
Eq. (ED). 

However, if SHDM particles are fermions like in the 
case of neutralinos, there is quite different effect which 
stops the core contraction at a much larger radius. This 
effect is the pressure of a degenerate Fermi gas. The 
maximum density of the core and hence the radius of 
the core can be derived from the equality of the Fermi 
momentum of a degenerate gas and the virial momentum 
of the constituent particles at the radius of the core r = 



V f = (37r 2 ) 1/3 (p c A«x) 1/3 - ™ X V C 



(25) 



where V c = v/ GM c /r c is the virial velocity at the bound- 
ary of the core, and M c = (47r/3)p c r 3 is the mass of the 
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core. For the density profile p(r) oc r~ 2 , the virial veloc- 
ity is the same at all r and we can take it for the whole 
clump with mass M and radius R. Using core radius 
from x c = (p/Spc) 1 / 2 we obtain 



-3/2 



(26) 



For a clump with mass M ~ 1 x 10 5 g, mean density 
p ~ 3 x 10 3 g/cm 3 and R ~ 3 cm, the radius of the core 
is x c ~ 1 x fCT 11 . 

In our calculations below for superheavy neutralinos, 
we will use for radius of the core Eq. (|26l) . 



C. Properties of superdense clumps from SHDM 
particles and numerical examples 

The detectability of the annihilation signal from DM 
clumps is determined by the following parameters: The 
mass m x of the superheavy neutralino and its mixing 
parameters, the density profile, e.g. an isothermal profile 
p(r) = p c (r I R c )~ 2 in case of a 'gravithermal catastro- 
phe', the dimensionless radius of the core x c = R c /R = 
{p/ipc) 1 / 2 and the maximal density p c . The important 
parameter, the mean density of a clump p, is found from 
the evolution of a primordial perturbation with initial 
amplitude Sh for a clump of given mass M . The mean 
densities are shown in Fig. [1] for different M and Sh- 
The maximum density, and respectively the radius of the 
core, are found as described above. 

As particular examples we consider three sets of clump 
parameters which we shall use in the next section for 
the calculation of the annihilation signal from superdense 
clumps. We discuss first an optimistic example for a bino 
with mass m x = 10 11 GeV. We consider clumps with 
M ~ 10 5 g formed from fluctuations with Sh — 0.09, 
marked by a star in Fig. [T] The density profile of such 
clumps is p(r) oc r~ 2 (we choose the isothermal value, 
which is close to the analytical and numerical results j3 ~ 
1.7 — 2) and the initial core radius x Cy i ~ 0.01. Such 
clumps have mean density p ~ 1.3 x 10 3 g cm -3 , radius 
R ~ 2.6 cm, and virial velocity v ~ 0.05 cm/s. For 
these parameters, the evolution time of the initial core is 
t e ~ 0.4^0- Fermi degeneracy mainly restricts the central 
density in this case, and the new core radius is R c ,2 — 
2.9 x lO" 11 cm (x c , 2 ~ 1.1 x lO" 11 ). 

The parameters of a more typical example (marked 
by a cross in Fig. P are M ~ 10 15 g and 5 H — 0.07. 
These clumps have the same density profile and initial 
core radius, p{r) oc r~ 2 and x Ct i ~ 0.01. But now p 
6.3 x 10~ n g cm" 3 



R ~ 1.6 X 10 s cm, v ~ 0.65 cm/s, 
and no singularity forms in the clumps. 

In the case of a higgsino with the same mass, m x — 
10 11 GeV, an optimistic choice of parameters corresponds 
to clumps with M = M d ~ 6 x 10 21 g. The density 
profile of such clumps is p{r) oc r~ 2 and the core radius 
x c .\ ~ 0.01. Also in this case no gravithermal catastrophe 




FIG. 1: The mean density p (in g cm -3 ) of DM clumps as 
function of the perturbation Sh in the radiation density on the 
horizon scale; solid lines from top to bottom are for for clump 
masses M = 1CT 10 , 10~ 5 , 1, 10 s , 10 10 , 10 15 , 10 20 , 10 25 , 10 30 , 
10 g. The dashed line is the bound on the clump density 
from primordial black holes overproduction with threshold 
S c — 0.7. The time of two-body gravitational relaxation inside 
the clump cores is less then the age of the Universe for clumps 
above the dotted lines for m x = 10 11 GeV. The star marks 
favourable parameters for annihilations, and the cross marks 
a typical example considered for comparison. 



develops. These clumps have p 
10 9 cm, and v ~ 600 cm/s. 



10 



-6 g cm 3^ 



R 



IV. ANNIHILATION 



A. Annihilation rate 



The annihilation rate JV ann 
clump is 



1 



4:Trr 2 dr n 2 (r)(a a 



of neutralinos in a single 



M 2 



8- 



i? 3 



where n(r) = p{r)/m x is the number density of neutrali- 
nos as function of the distance to the core of the cloud. 
The function 5* was determined in Ref. [2l| and depends 
on the distribution of DM in the clump. In particular, 
the function is S = 1 for the simplest case of an uniform 
clump and S ~ 4/(9x c ) for an isothermal profile p oc r~ 2 
with a small core size, i c « 1. 

The resulting flux Ii of particles of type i = N,j,u 
from DM annihilations summed over all DM clumps in 
the Galactic halo is given by 



1 



1 dNi 

7i v dx 



(28) 



where dNi/dx is the differential number of particles of 
type i produced per annihilation with energy E = xm t . 
We calculate these spectra as described in Ref. @, [29| 
for the case of a non-supersymmetric evolution of the 
fragmentation functions dNi/dE. 
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The function H contains the information about the 
smooth DM distribution in the halo, 



d( sin ( 



A \i^Mfl, (29) 



where £ is the fraction of DM in form of neutralino 
clumps, n c \(r) = ^p^/M is the number density of clumps 
at distance s from the Sun along the line-of-sight (l.o.s.), 
and £ is the angle between the direction in the sky 
and the galactic center (GC). Finally, r max = (R^ — 
Tq sin 2 C) 1 ^ 2 + r & cos£ is the distance to the border of 
the DM halo of radius and r = 8.5 kpc is the dis- 
tance of the Sun to the Galactic center. 

As distributions of the DM in the galactic halo we use 
the Navarro- Frenk- White profile [32[ , 



Ph(R) = 



(R/R s )«(l + R/R s )< 



(30) 



with a = 1, scale radius R s — 20 kpc, Rh = 200 kpc as 
the size of the DM halo and Ph(^o) = 0.3GeV/cm 3 as 
the DM density at the position of the Sun. 

Annihilations may proceed in the clumps with differ- 
ent parameters because the a priori unknown position 
and height of the putative spike in the spectrum of per- 
turbation. The formation and the properties of super- 
dense clumps were considered in Paper I [15j . Following 
the formalism of Ref . [l5| , the clumps density follows as 
shown in Fig. [T] For illustration we use in our calcula- 
tions the three sets of clumps parameters presented in 
the previous section. 



B. Cross section 

Analytical approximations for the annihilation cross 
section of a neutralino valid in the limit A/susy 3> mz 
were presented in Ref. jTTj using lowest order pertur- 
bation theory. We use as annihilation cross section for 
both binos and higgsino the subprocesses XX ~^ ZZ and 
XX ~^ W + W~ in the case of a higgsino, 



7)~2x 10~ 42 m 



li 



2 cm 3 s~ 



(31) 



Annihilations into Zh° and ZA can increase the cross- 
section, depending on the values of the SUSY breaking 
parameters. 

Since the relative velocities of SHLSPs in DM clumps 
arc small, /J = u <IC 1, factors g 2 / (3 or ln(g 2 / '/?) can 
lead to a break-down of perturbation theory. This effect, 
first studied by Sommerfeld for Coulomb interactions, 
was generalized in Ref. 30] to the exchange of massive 
non-abelian gauge bosons relevant for neutralino annihi- 
lations. In the case of a wino or higgsino, the small mass 
splitting 5m — m x — m x ± between the lightest neutralino 
and the lightest chargino means that the charginos pro- 
duced in XX xtxi have the same small velocity as 
the neutralinos. Therefore multiple photon, Z and W ± 
exchange between the charginos becomes important. 




E[eV] 



FIG. 2: The maximal fluxes Ii(E) of photons, nucleons and 
neutrinos from neutralino annihilations in Galactic halo to- 
gether with experimental data for a neutralino with 10 11 GeV. 



The resulting enhancement of the annihilation cross 
section can be calculated non-relativistically and is, ne- 
glecting bound-state effects, characterized by two param- 
eters [3l[ : The ratio e — mw jm x determines, if the an- 
nihilation proceeds in the Coulomb (e <C 1) or in the 
Yukawa (e 1) regime, while the ratio x = g 2 s /f3 of 
the squared effective coupling constant and the veloc- 
ity determines, if factors g 2 g/f3 or ln(g 2 s /f3) lead to a 
break-down of perturbation theory. Here, the effective 
coupling constant <? ff includes all pre-factors in front of 
the Yukawa potential, as e.g. mixing factors. 

The Sommerfeld factor 1Z as ratio of the perturbative 
and non-perturbative annihilation cross section is given 
in the Coulomb case by 



n 



o-port 1 - exp(-r?) 



(32) 



with -q = ±(? 2 ff /(2/3). Using in the case of a higgsino the 
parameters for the optimistic example given in Sec. lIV K\ 
we find a strong enhancement of the annihilation rate, 
K - 10 10 . On the other hand, hig gsino are relatively 
tightly coupled to the thermal plasma, leading to a rela- 
tively large value of Md- Thus there is no gravithermal 
catastrophe for higgsino clump and we use 5=1. 

For a bino, the Sommerfeld effect is not effective, 
TZ = 1. Using the above parameters for the optimistic 
example, we find S ~ i/(9x c 2) ~4x 10 10 . Hence, anni- 
hilations of binos in clumps formed during the RD epoch 
can be enhanced by the factor S ~ 4 x 10 10 compared 
to the annihilation signal computed for a smooth DM 
distribution inside clumps. 

The maximal fluxes of photons, nucleons and neutrinos 
from a SHLSP with m x = 10 11 GeV allowed by cosmic 
ray data are shown in Fig.[2]togcthcr with upper limits for 
the neutrino fluxes from the Pierre Auger experiment [33[ 
and RICE as well as the expected 90% sensitivity after 
five years data-taking of ICECUBE [34| . For the values 
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given in the optimistic example for a bino, the flux shown 
in Fig. [5] has been rescaled by the factor 10 -14 , while 
for a higgsino the flux was rescaled by 10 -5 . Thus the 
fluxes of secondaries overshoot the experimental data by 
many orders of magnitude for the most optimistic scenar- 
ios. This offers the potential to test different DM masses 
as well as different scenarios for the formation of super- 
dense clumps. At present, annihilations of SHLSPS are 
mainly restricted by experimental limits on the photon 
fraction 35| and a galactic anisotropy of the UHECR 
flux (36j . However, in the future neutrino searches at 
lower energies by a km 3 neutrino telescope as ICECUBE 
may become competitive. 



V. CONCLUSIONS 

We have studied the properties of SHDM clumps in two 
different cosmological scenarios, one with power-law and 
one with spiky density perturbations. As superheavy DM 
particles we have considered the superheavy neutralino, 
which properties were studied within superheavy super- 
symmetry in Ref. [TTj . 

For standard power-law fluctuations, SHDM clumps 
are formed in the DM dominated epoch in hierarchical 
structures, when a small clump is hosted by a bigger 
clump, which in turn is submerged into an even bigger 
one, etc. Small clumps are tidally disrupted in such struc- 
tures and only a small fraction of the clumps survives. 
The surviving clumps can be further destroyed by tidal 
interaction in the Galaxy. 

In contrast, clumps in spiky density perturbations are 
born in the RD dominated epoch, when hierarchical 
structures are not yet formed. They evolve as the single 



isolated objects without tidal destruction and merging. 

For the standard cosmological scenario the density of 
SHDM particles in a clump is low and the annihilation 
signal is weak, caused by the too small annihilation cross 
section. These clumps can be detected only when clumps 
are passing by gravitational wave detectors. 

In superdense clumps, an isothermal profile may con- 
tinue up to very small radius of the new core, if the clump 
went through the "gravithermal catastrophe" . The den- 
sity of particles in the core and nearby is very large and 
this enhances the annihilation signal. Another reason 
for the increase of the annihilation signal in comparison 
with clumps of EW mass particles is the small M m ; n in 
the mass distribution of the clumps. The cutoff M m i n is 
smaller for a bino, but a higgsino gains from the Sommer- 
fcld factor, which increases the annihilation cross section. 

As a result, we found that the annihilation rate of sta- 
ble superheavy neutralinos may be large enough to be 
detectable, if primordial density perturbations are spiky. 
Hence the search for photons or a galactic anisotropy in 
UHECRs 36j as well as the search for UHE neutrinos of- 
fers not only the potential to identify the DM candidate 
but also to learn about the inflationary potential. 
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